r k 



Specify the (i) chiller make and model, (ii) the frequency range of the spectra = [F range lo , 
F range hi ] Hz, (iii) electrical line frequency, (iv) importance levels 



Get all the vibration spectra for the chiller from the historical database 



Calculate the motor frequency, Mjj , for every spectrum Sj and scale the 
frequency axis of Sj by its motor frequency, 



<n 



Find the spectrum with the smallest no. of frequency lines = m 



<3 



N bw =1.5x(F rangehi -F rangel0 )/(mxM/) <jT| 



For i 



= 0 to 2m 



fi = low + « N bw /3)/M 7 (T 



Collect vibration data from all spectra in the frequency range [f { -2N bwf f i +2N bw ] 



Perform outlier detection and removal using Caroni and Prescott's method 



Calculate kernel density p(x \F=f t ) and its cumulative probability P(x \F=f { ) ^10 



Calculate the value of x = such that P(Xp \F=f i )=p- / 
Xp is the 100j8% amplitude limit at frequency /j. 



Figure 1 




FIGURE 2 
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Figure 3 



Given dataset 









Set «out = 0 




1 





21 



> For/= l,2,...,rdo 



Compute the average of elements for the dataset X 



Compute the matrix of sums of squares of outer 
products, A p for the dataset X 




Find i* extreme observation . in dataset X 

e,t 



Compute the I th Wilk's statistic D. for multivariate data 



Compute the / th critical value y. 



27 




Set n = j 

out 



Remove the i* extreme observation 
x . from the dataset X 



30 



( Outliers: {x„ „ x .... x } 

V *■ e,V e,2' ' e flout > 
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Figure 4 



0.03 




0.47 0.48 0.49 0.5 0.51 0.52 0.53 

fx motor frequency 



Figure 5 
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Kernel name Definition 



Plot 



1 Gaussian 



2. Epanechnikov k(u) = |(1 - u 2 ) l{| M |<i) 



3. Biweight k(u) = (1 - u 2 ) 2 1, M<1 , M 



4. Triangle k(u) = (1 - |w|) l {M< i) M 




Note: l(| U |<u = 



|l ifM^l 
0 otherwise 



FIGURE 6 
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FIGURE 7 




FIGURE 8 
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FIGURE 10 
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FIGURE 11 




FIGURE 12 



